THE FREQUENCY DETUNING CORRECTION AND 
THE ASYMMETRY OF LINE SHAPES: 



THE FAR WINGS OF H 2 0 - H 2 0 


Q. Ma 

Department of Applied Physics, Columbia University 
and Institute for Space Studies, Goddard Space Flight Center 

2880 Broadway 
New York, New York 10025 


and 

R. H. Tipping 

Department of Physics & Astronomy, University of Alabama 
Tuscaloosa, AL 35487-0324 


/ 


2 


ABSTRACT 

A far— wing line shape theory which satisfies the detailed balance principle is applied 
to the H 2 0 — H 2 0 system. Within this formalism, two line shapes are introduced, 
corresponding to band— averages over the positive and negative resonance lines, 
respectively. Using the coordinate representation, the two line shapes can be obtained by 
evaluating 11— dimensional integrations whose integrands are a product of two factors. One 
depends on the interaction between the two molecules and is easy to evaluate. The other 
contains the density matrix of the system and is expressed as a product of two 
3-dimensional distributions associated with the density matrices of the absorber and the 
perturber molecule, respectively. If most of populated states are included in the averaging 
process, to obtain these distributions requires extensive computer CPU time, but only have 
to be computed once for a given temperature. The 11-dimensional integrations are 
evaluate using the Monte Carlo method, and in order to reduce the variance, the 
integration variables are chosen such that the sensitivity of the integrands on them is 
clearly distinguished. Numerical tests show that by taking into account about 10 7 random 
selections, one is able to obtained converged results. We find that it is necessary to 
consider frequency detuning, because this makes significant and opposite contributions in 
the two band-averaging processes and causes the lines to be asymmetric. Otherwise, the 
two line shapes become symmetric, are the same, and equal to the mean of the two shapes 
obtained including the frequency effects. For the pure rotational band, we find that the 
magnitude of the line shape obtained from the positive line average is larger than that 
obtained from the negative fine average for u > 0 and vice versa for u < 0, and their gap 
increases as the frequency displacement from the line center increases. By adopting a 
realistic potential model and optimizing its parameters, one is able to obtain these two line 
shapes and calculate the corresponding absorption coefficients which are in good agreement 
with laboratory data. Also, this same potential yields good theoretical values for other 
physical properties of the dilute H 2 0 gas. 
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I. INTRODUCTION 

In previous papers, we have presented the theory for the calculation of far— wing line 
shapes and the corresponding absorption coefficients for interacting pairs of molecules. 1 ' 3 
Assuming only the binary collision and quasistatic approximations, we have shown that by 
using the coordinate representation to describe the orientation of each molecule before and 
after the transition, we are able to reduce the problem to the calculation of 
multi— dimensional integrals. The dimensionality of the integrals depends on the type of 
molecules involved; specifically for two linear molecules (e.g. C0 2 — C0 2 or C0 2 — N 2 ) the 
dimensionality is 7, while for one asymmetric top and a linear partner (e.g. H 2 0— N 2 ) it is 
9. For the first case we were able to obtain converged results with a sophisticated 
interaction potential using conventional integration methods. 1 However, for latter case, we 
had to use the Monte Carlo method. 3 This system is important in atmospheric 
applications, where we have shown that the absorption (known as the "foreign continuum") 
is in good agreement with experimental results. 3 Because the absorption by H 2 0— H 2 0 
pairs (the "self continuum") is more important, we would like to extend our theory to the 
case of two interacting asymmetric tops; the dimensionality in this case is 11, thus 
implying a big challenge to obtain converged results. In a previous paper, we have shown 
that by considering an interaction potential containing cyclic coordinates, the 
dimensionality is reduced to 7, and one can obtain converged results. 2 

One goal of the present paper is to remove this restriction and calculate the 
far— wing line shape for H 2 0— H 2 0 using the most general interaction potential. To 
accomplish this goal, it is necessary to modify the Monte Carlo routine used previously 
such that the sensitivity of the integrand on the integration variables is clearly 
characterized. A second goal is to investigate the asymmetry of spectral lines; in order to 
do this, one has to carry out band averages in a more sophisticated way and consider the 
frequency detuning effect in the line shape calculations. Based on the present work, one 
can conclude that the band— average line shapes are asymmetric and, in addition, one would 


r 
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expect different line shapes for different bands. 

In order to reduce any unphysical effects in calculated results, we carry out all 
numerical calculations based on formulas which satisfy the detailed balance principle 
exactly and have a higher accuracy in the short-time limit. 4 As expected, this increases 
the difficulty because the formulas become more complicated. Thanks to great advances in 
computers made in recent years, we are able to overcome this obstacle successfully 
Finally, we utilize a realistic interaction potential applicable for the dilute H 2 0 gas that 
not only gives good agreement between theory and experimental temperature-dependent 
second virial coefficients and differential scattering cross-sections, but also yield absorption 
coefficients that are in good agreement with the experimental self continuum data. We 
note that because the theory has a sound physical basis, we expect that theoretical results 
calculated at lower temperature, which would be extremely difficult to obtain 
experimentally because of the lower vapor pressure attainable, are valid. 

The paper is organized in the following way. In Sec. II. A, we present the 
expressions for the correlation functions and the spectral densities which satisfy detailed 
balance exactly. In Sec. II.B, we introduce two band averages and the corresponding line 
shape functions. The application of the coordinate representation for the system consisting 
of two asymmetric tops is presented in Sec. II. C. The main purpose is to introduce the 
density matrices of the absorber and the perturber molecule, which are three dimensional 
distributions, and to show how to calculate them. This is the most costly calculation in 
the present study and the results obtained are used to get the line shapes later. In Sec. 

II.D, we discuss how to calculate the frequency detuning correction. The method is similar 
to that used to calculate the line shapes without this correction, except one has to 
developed a new technique to deal with the integrand containing a derivative. The 
necessary generalization of the Monte Carlo method to 11-dimensional integrations is 
discussed in Sec. II. E. In Sec. II. F, we discuss a realistic interaction potential which gives 
good agreement with various molecular data measured in the dilute gas phase, including 
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the absorption coefficient. Then in Sec. II. G, we apply the theory to obtain two line shape 
functions numerically for several temperatures and discuss some of their general features. 
We use these results to calculate the corresponding absorption coefficients and compare 
those obtained in the 300 — 1100 cn" 1 spectral region for the room temperature with 
experimental results. In Sec. Ill, we discuss the results obtained and the conclusions drawn 
from the present study. 

The present work, together with the previous calculations for simpler systems, 
constitute a general formalism in which one can obtain from first principles the far— wing 
line shapes for any colliding pair for which a realistic potential is available. Conversely, if 
experimental data for the line shape or the corresponding absorption coefficients are 
available, one can use this formalism in order to test the accuracy of the potentials. 

II. THE GENERAL FORMALISM 

A. Symmetric Correlation Function and Spectral Density Which Satisfy the Detailed 
Balance Principle 

The absorption coefficient a(u > ) of a gaseous sample with a unit volume is given by 

o_2 

where n a is the number density of absorber molecules and the spectral density, F(w), is the 
Fourier transform of the correlation function C(t). One separates the total Hamiltonian H 
into two parts: one commutes with the internal coordinates of the molecules while the 
second does not. We note that this distinction of H usually coincides with the division of 
the interaction into two parts: and V ani , the isotropic and the anisotropic interactions, 

respectively. Accordingly, the total Hamiltonian H is decomposed as 

H = H a + H b + + V ani = H 0 + V Mi , (2) 

where H a and H b are the unperturbed Hamiltonians of the absorber and the bath 
molecules, respectively. In practice, for atmospheric applications where the gas pressures 
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are low, one can introduce the binary collision approximation and focus on a much simpler 
system consisting of one absorber and one bath molecule. For simplicity, we don’t 
introduce new symbols for the two-molecule system. 

It has been shown that one is able to express C(t) in the t -* 0 limit, which is valid 
to order t 2 , as 4 


C(t) = n b v Tr{(e^ at v^ p) f P^e ^ani [ e 1 W( e “ ^ (3) 

where n b is the number density of the bath molecules, v = Tr[e"“^ H °]/Tr[e _ ^ H ] and the 
Liouville operators L & and Z< ani corresponding to H a and V ani , respectively, have been 
introduced. For later convenience, we can introduce a symmetric correlation function C(t) 
(= C(t + iti/3/2)) defined by 


C(t) = n b <-Tr{(ei i “*v^4/‘^ t 

« Pis o e_ ^4pI pj P pj)] e - (4) 

From these expressions, it is easy to verify that C(-t) = C(t + ih$) and C(-t) = C(t) 
which guarantees both of them satisfy the detailed balance principle. 4 ’ 5 

By expressing C(t) explicitly as summations over indices i, j, i’ and j’ where each 
represents all the quantum numbers necessary to specify the energy levels of the absorber 
molecule, one is able to write C(t) as 


C(t)= l n b „ Tr b {<j|pf /i +4^|i> 

ij i’j’ 

<* I Piso e ~~^^ &ai e“ lVanit |i’> <i’| p\ MPllj’> <j’|e iV * nit e-^ani) j>}. (5) 


In the above expression, = Ej — Ej and the the subscript b of trace denotes the trace 


over the remaining variables including all magnetic quantum numbers. We assume that 
the interaction between two molecules does not depend on their vibrational quantum 
numbers. As a result, in Eq. (5) the vibrational quantum numbers of i and i’ are identical 
and those of j and j’ are also identical, but the former could differ from the latter. By 
choosing the Z axis of the space— fixed frame along the separation between the two 
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molecules, the interaction V ani (r, fl a , fl b ) depends on the orientations of the two molecules 
represented by fl a and ft b , respectively, and on r, the distance between the centers of mass, 
which can be considered as a parameter since the translational motion is treated classically. 
In Hilbert space associated with the internal degrees, one denotes the eigenvalues and 
eigenvectors of V Mi (r, fl a , fl b ) by G c and | <>, respectively; thus 

^ani( r > ^a> ^b)l = G^(r) | (6) 

Then, by performing the Fourier transform of the correlation function C(t) and carrying 
out an integration over r (i.e., the classical ensemble average over r which is valid within 
the quasistatic approximation), one is able to obtain the symmetric spectral density F(u), 

*■(") = ; I I ( 7 ) 

. . .) . ) 

1J 1 J 

The explicit expression for Jy.pj^w) is given by 

= l H Ct) M l l l <T7| ji b > <ji b l4AlV^|ii b > <ii b IC> 

^ { m } i b i b 


* <C|i’i b > <i’i b |V^PlMPllj’i b > <3%\v>, (8) 

where i b and i£ are indices used to specify the states of bath molecule and {m} indicates 
the summation over all magnetic quantum numbers. In the above expression, H^(w) is 
defined by 


h ;t,( w ) = n b * 47r2f c 




g — /?Viso(rc) — /?[ G <;(rc)+ G r 1 ( r c)]/2 


(9) 


where G^(r) denotes jp[G;(r) - G T) (r)], and r c are roots of the equation 

G ;( r c) “ G ti( r c) = u - (10) 

We note that in Eq. (7), the summation indices i, j, i’ and j’ exclude their magnetic 
quantum numbers since the summation over them has been carried out. The functions 
*ij.i,j,(w) introduced above are symmetric for the exchange of indices {i j} « {i’j’}. 
Meanwhile, with respect the exchanges {i j} -* (j i) and {i’j’} -» {j’i’}, it is easy to verify 
that — w) = Jij.i,j,(w). 4 In terms of F(oi), the absorption coefficient a(w) becomes 
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0—2 

= gjg n a w sinh(/3?ia>/2)F(w). ( 11 ) 

Finally, it is worth mentioning that F(w) = e - ^"/ 2 F(a/) and F(-w) = F(w). The latter is 
required by the detailed balance principle. 

B. The Averaged Line Shapes with the Frequency Detuning Correction 

We consider a band consisting of transitions between states with two vibrational 
quantum numbers specified. For simplicity, we designate a pair of i and j (i.e., a line) by 
the symbol k. If necessary, we use symbols k + and k. to distinguish a positive resonance 
line with Ej — E A > 0 and a negative resonance one with Ej — E i < 0. Then Eq. (7) can be 
expressed as 

F M = \ {J \ * k+; J w “ + w n)l + W * k _; n [ w + ~ w n)]}> ( 12 ) 

k+ n k- n 

where (= | u^l ) is positive. We note that in Eq. (12) we have not explicitly 
characterized the symbol n by + or — yet because the band of interest has not been 
specified. For the pure rotational band resulting from transitions without changing the 
vibrational quantum number, both the symbols n in two terms of Eq. (12) could be n + and 
n. because there is no link between the choices of k and n. Meanwhile, for a vibrational 
band, k + or k. is always associated with n + or n_, respectively. In this case, i and i’, and 
also j and j’, share the common vibrational quantum numbers and, in general, these two 
different vibrational quantum numbers determine whether k (i.e., a pair of i and j) and n 
(i.e., a pair of i’ and j’) belong to the positive or the negative resonance. Therefore, for the 
vibrational bands the symbol n in the first term of Eq. (12) should be understood as n + and 
those in the second term as n_. This custom is applicable for the equations following. 

On the other hand, by introducing two common line shape functions J + (u;) and X.(u) 
one can express F(w) as 

= x ^ ^klflkl 2 1 u k + ) + + 1 w k+ )2 ~ X i u + w k*)]}. ( 13 ) 

k* 


/ 


9 


where (= py} is the reduced dipole matrix element and p k = In order to find 

expressions for <T + (o/) and X.{u), one performs the Fourier transform for Eqs. (12) and (13) 
and compares the results to their time-domain versions. As a result, the following 
equation has to be satisfied: 

k+ n k- n 

= ^A(“)]^kl^l 2 ^*‘+ (14) 

k+ k- 

where for simplifying the notations we use the symbols x . (t), x . (t), &[K J + (u;)], and 

k+jn k-»n ^ 

^[^5 *-(<*>)] to represent the Fourier transforms of x (a/), x (u>), \ X+(v), and 

k +) n k-jn ^ ^ 

X.(uj), respectively. Furthermore, one assumes that one can separate Eq. (14) into two 
equations, one associated with a summation over the positive resonance lines and the other 
over the negative resonance lines, 


■’TsjAMI-X Rj 1 ) e^* +W ” )l / I 


k + n 


(15) 


s + 


and 


k- n s ♦ 

Then, with Eq. (15) it is easy to obtain an expression for /T + (u))] which is valid in the 
short-time limit: 

>.(«)] l l-W e‘* Wk * + “» ) " “ ■*, (17) 

k+ n 

where N is the normalization factor defined by 

N = I />kl/*kl 2 > ( 18 ) 

k + 

and Z) is the average positive resonance frequency defined by 

w = lSrX ^kl^kl 2 Vk*- (I 9 ) 

k + 

By performing the inverse Fourier transformation, one is able to derive the expression for 
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X+(u) valid at | cj\ — ► oo 


x*( u ) = tJ2 bl I - + u a) + u ]• 

k+ n 

(20) 

Similarly, one can derive 


Xlu) = U* J \ “ ^n) “ u )• 

k- n 

(21) 


In Eqs. (20) and (21) the arguments of functions depend on the summation indices 
and this results in difficulty to obtain X+(u) and X _(w) directly. In practice, one prefers to 
derive an expressions for X + (u>) and X.(u>) inwhich the summations are performed over 


functions whose arguments are independent of the summation indices. For this purpose, 
the frequency detuning approximation must be introduced. As an example, we consider 
Eq. (20). With the Taylor series expansion of x . [u> - + u a ) + u> ] over u, one can 

k*>n 

approximate X.(u>) as 

X + (u>) = X+(u) + u) 2 ^ ^ J [ui - |(a^ + + w n )] x! . (^), ( 22 ) 

^ k+jn 

k + n 


where 


Euj7 ^1 


k* n 


(23) 


and x! . (w) = dy . ( u)/du ;. In the left side of Eq. (22), the first term XJoj) comes from 

k+>n k+>n v ' 


simply ignoring the frequency detuning and the second term is a correction. We note that 
in our previous study, 6 a shift parameter was introduced to treat the effects from the 
frequency detuning. However, it is better to calculate the contributions from the second 
term directly because it turns out that, except for some simple cases, the previous method 
could introduce numerical errors. Similarly, one can approximate ?.(w) as 

X^v) = X.(u) + W 2 J}- J J [- w + - w n )] ^ j n ( w )> (24) 

k- n 

where 

i-(“) = 1 

k- n 


( 25 ) 
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Using the symmetries y . (-w) = x, . (w) and x . (-w) = X . M (the latter is 

k-jn- k+in+ k->n + k+>n-. 

appropriate for the pure rotational bands only) mentioned above, one can show that 

iji-u) = i + (w), (26) 

and 

Xl-d) = > + (u»). (27) 

If one does not distinguish the two shapes X+( u) and X.{u) in the expression for F(w) 

given by Eq. (13) and replaces them by only one shape X(u>), one can pursue a similar 

derivation for X(u) and obtain 

X(u>) = u 2 H y k . n [u; - + w n )] /\p k | /^ | 2 . (28) 

1 n k 

A simple version X(u) defined by 

X(u) = J 2 H x (u) / l |^kl 2 , (29) 

k n k 

can also be introduced. In this case, both X(u) and X(u)) become symmetric and, in 
addition, up to the first-order approximation, X(u) is the same as X(u) because there is no 
net contribution to X(u) from the first derivative term of the Taylor series expansion of 

x k]n W-^ + » a )]. 

In summary, we note that within this formalism no matter which functions are 
chosen for the line shape, in terms of them the symmetric spectral density F(w) always 
satisfies the detailed balance principle. The formalism outlined above was developed 
several years ago.® However, at that time, except for the simplest system, such as C0 2 — 
Ar, to calculate converged line shapes for systems of interest in atmospheric applications 
was formidable. The main obstacle was the calculations involving a diagonalization 
procedure of the anisotropic potential which exhausts computer resources very quickly. 

C. The Coordinate Representation 

Recently, we have developed a formalism based on the coordinate representation in 
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which the eigenfunctions of the orientations of the system are chosen as the complete set in 
Hilbert space. 1 ' 3 The advantage of introducing this representation is that the 
diagonalization of the potential becomes unnecessary and the main computational task is 
transformed to the carrying out of multidimensional integrations. For systems consisting 
of two linear molecules, or one linear and one asymmetric top molecule, or two asymmetric 
top molecules, the dimensionality is 7, 9, and 11, respectively. In addition, we have shown 
that using the Monte Carlo method, one is able to evaluate up to the 9— dimensional 
integrations required for systems such as H 2 0-N 2 . 3 “Combined with techniques developed 
recently to handle sophisticated potential models, 2 one is able to implement realistic 
potentials for these systems and derive accurate, converged results for the far— wing line 
shapes and the corresponding absorption coefficients. 

In the present study, we are interested in a system consisting of two H 2 0 molecules. 
In order to reduce any unphysical effects, we base our study on formulas which satisfy the 
detailed balance principle exactly and have a higher accuracy in the short-time li mi t 
Besides extending the Monte Carlo method to evaluate 11-dimensional integrations 
required for line shape calculations, we would like to answer some questions in depth. We 
want to know whether the line shapes are asymmetric and if they are, to find out the origin 
of the asymmetry. In other words, we want to know whether it is necessary to introduce 
two line shapes instead one, and why they differ from each other. As expected, in 
comparison with our previous studies, this introduces extra difficulties. Fortunately, 
thanks to the coordinate representation and the Monte Carlo method as two powerful tools 
to perform averages, and also to great advances in computers made in recent years, we are 
able to overcome these obstacles and we make significant progresses. 

The details about the coordinate representation and the Monte Carlo method have 
been presented previously 1 ' 3 and are not repeated here. We only report new features. 

Since we want to carry out band averages in a more sophisticated way, we have to 
introduce a positive and a negative resonance dipole operator in the Hilbert space of the 


13 


absorber molecule denoted by and /z^J, respectively. The former is defined by 

I <»lMmlj> |i><j|, (30) 

Ej > Ei 

and the latter by 

= l <i|0jj> |ixj|. (31) 

E j< E i 

We note that for the vibrational bands j denotes a higher vibrational quantum n um ber 
than i in Eq. (30) and vice versa in Eq. (31). Then, with Eqs. (8) and (23), one can rewrite 
J + (w) as a summation over £ and rj 

X^u) = a/ 2 ^ J G ( * Ct) , , (32) 

C*i 

where G*^, are defined by 

G Kt|) = \ {<(WPbPllxZ Pi \v>)*<(W(^Pltt m pl\v>, (33) 

m 

and defined by 

G(C 1 ,> Pi \v>)*<(\Jpb f>t P* pt\v>> (34) 

m 

for the pure rotational bands and for the vibrational bands, respectively. For simplicity, 
we will only present formulas applicable for the pure rotational bands and simply mention 
differences between them and their vibrational analogs. In the expression for X+(uj) given 
by Eq. (32), the summation terms are products of H^(w) and G^, . The former are 
functions of u and their values depend on the interaction potential between the two 
molecules. The latter are common for all frequencies and their values are independent of 
the potential. With the coordinate representation, no matter how complicated the 
potential is, to calculate values of H^(w) is straightforward since the potential is a 
diagonal operator. On the other hand, to obtain the G^, involves a lot of calculations, 
because they contain the density matrices which are differential operators. Fortunately, for 
a given temperature, these calculations need to be done only once since results obtained are 


/ 
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applicable for all potential models. Thus, we can calculate them and to store them in an 
input file. 

In the coordinate representation, the basis functions | £> are nothing but the direct 
product | £(ft a — n a( -)> ® | £(n b — n b; )> where the notations of | £(O a — fl a( .)> and | 5(ft b — 
fi b; )> are used to represent specified orientations of the absorber and the perturber 
molecules, respectively. Accordingly, one can separate the dependence of on the 
absorber and on the perturber molecules and express it as the product of G a( ^ } and 

Cn) > 

g k„, = I «*(«. - «WI/>1 4l «(0. - n ar| )>)*<«(n s - n^j \,\ ^ „fl «(n, - aj> 

m 

« |<«(n b -o b j|v^|tf(n b -n ta| )>|> 

- * ®b(Cn) ' ( 35 ) 

For the vibrational bands, the dipole moment operator /z m in Eq. (35) is replaced by 

In comparison with G a( ^ , the expressions for G b( ^ , are simpler because they do 
not contain the dipole operator. For the linear, symmetric top, and asymmetric top 
molecules, the explicit expressions for G b( ^ have been presented and the corresponding 
profiles have been discussed in our previous papers 1 " 3 and we do not repeat them. With 
respect to G a( ^ , one has to derive the corresponding expressions valid for the linear, 
symmetric top, and asymmetric top molecules, respectively. We do not present all of 
them, rather only the last and the most complicated one applicable for H 2 0. It is well 
known 7 ’ 8 that the wavefunctions of H 2 0, | jrm>, can be expressed in terms of an expansion 
of symmetric— top wavefunctions | jkm>, 

|jrm> =^u£ r |jkm> 

= l U kr K- 1 ) m _ k A 7)1. (36) 

k 

where ^(a, /3, 7) (= e - * 1110 ^ is the rotational matrix. With Eq. (36), one 

is able to express G*^, as, 
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II ^KK’ ^K,K^ a < Cn> » A Cn> 1 Cn> 


(37) 


L KK’ 


where a ( > A ^j) > 7 ( ;n) 416 three Euler angles used to represent a rotation 

resulting from two successive rotations, i.e. , 

^■^(Cn)’ Ac*!)* ^(Cn)) “ ^ ( a C’ Ap 7^) ^(^n* Arp 7q)> (38) 

and the summation index L = 0, 1, 2, • • • ; both indices K and K’ run from — L to L. In 

the above expression, the coefficients * are given by 


A L+ — * ( i\l+I<+K 

a kk> - ^ 


I I (2j 1 +l)(2j J +l)vI7S7 e- WBa..^+E0»r J] 

1 2 


j l r 1 J2 t 2 


x I ^ (- l) j r +j 2 ^jl+l)(2j>+l) e M E (j’^i)+ E (j 2 .n)] W(j’ j 2 j’ j^lL) 

{ii r l}+ J2 r 2 

* (I (- !) k KV U k-K r, C (Ji h L, k K-k K)} 

k 12 

. > > 

« {£ (- l) k ’ui,' +K , r , uj? r> C(ji j; L, k'+K- -k’ K')} 

k ’ , ' ' 

* <1 U n|r, U n|rl C 0 1 «■ 0 « <1 U n’r, U n|r S % 1 & 0 “>»• ( 39 > 

“l n 2 

where E(j,r) are the energies of the state labeled by the quantum numbers j and r, g^ is its 
nuclear spin degeneracy factor, C(j t j 2 L, k K-k K) is a Clebsch-Gordan coefficient, W(jJ 
j 2 Ji L) is a Racah coefficient, and the summation over jj and t\ indicated by a symbol 
{j’i T iK is limited to a range with EQ’^i) > EQ^). We note that for the vibrational bands, 

T + 

the expression for A KK , is similar to that given above except that the states labeled by j[r\ 
and j 2 r 2 and the states labeled by jjTj and j 2 r 2 belong to a higher and a lower vibrational 
level, respectively. In addition, another limitation of E(jJrJ) > E(j 2 r 2 ) is also enforced in 
the summations. 


With Eqs. (37) and (38), one can conclude that G* ( ^ T)) are three-dimensional 
distributions over three Euler angles used to represent a rotation of the molecule from the 
initial orientation to the final one, labeled by £ and 77 , respectively. With respect to ^.(w), 
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one can introduce G^t c*i, an( * wr * te down similar equations to Eqs. (32) and (35). 
Meanwhile, G a( ;t)) can be given in terms of a£ k , and D^ K ,(a ( ;t)) , /?, ;t)) , 7( ;t)) ) in the 
same way as shown by Eq. (37). The expression for Aj , is almost identical to A*;* shown 
by Eq. (39), except that the summation over j’j and t\ is limited to a range with E(j’ 1 r’ 1 ) < 

Although the calculations of A^ K , and A^ K , are straightforward, there are many 
summation loops involved. Usually one introduces a cut-off to exclude less populated 
states. It turns out that as j max increases, not only more CPU time is required to calculate 
each A kk , and A KK , because the ranges of loops become larger, but also the number to be 
evaluated increases quickly. Fortunately, one does not need to calculate all of them since 
some are identical and others are zero. For H 2 0, due to the symmetry of Uj" , all the 
coefficients A KK , and A KK , are zero unless their indices K and K’ have the same evenness or 
oddness. In addition, for the non-zero coefficients there are symmetries A^ + , = A L+ ,, 

KK K— K 

^kk’ = ^-k-k” anc * ^KK’ = ^K’K‘ a result, if one introduces a cut-off = 23 (which 
is the highest angular quantum number of the initial states listed in the pure rotational 
band of the HITRAN 92 database 9 ), there are 18424 values of A*L, and 18424 of A[L , 
needed to be evaluated. If one uses an even higher cut-off j max = 26, these numbers 
become 26235. By utilizing a dozen workstations, we are able to manage the latter in less 
than two days. We note that to obtain these coefficients is the most costly calculational 
part in the present study. 

After all A KK , and A KK , are available, we can easily calculate G a( ^ and G a( ^ , 
which are three-dimensional distributions over the Euler angles a ( ^ , /? ( ^ , and 7( ( - T)) . 

In cases where no confusion results, the subscripts of these Euler angles are omitted. 
However, it is better to express them as distributions over their two sensitive variables J3 
and u (= (a + 7 )/2) and one insensitive one v (= (a - 7 )/2). The explicit expression for 
G a( used in the numerical calculations is given by 


! 
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G;. w (A«. T )-n a k«. «„■ 

L KK’ 

x {cos[(k+k’)u] cos[(k— k’)v] - sin[(K+K’)u] sin[(K- e’)v]}, (40) 

where the ranges of the indices L, K, and K’ are from 0 to 2j mgy , from 0 to L, and from - L 
to L, respectively; e KK , = 1 for K = 0, and e KR , = 2 otherwise. The expression for 
is similar to Eq. (40), except a replacement of A^*, by a£ k ,. However, due to the 
symmetry A RK , = A K , £ mentioned above, G* ( ^ (0, u, v) does not differ from G* ( ^ , (0, u, 
v) significantly. In fact, it is easy to show that 

G a( cn) (0» u > v ) = G a( cip (0. u > “ v )- ( 41 ) 

This means that with respect to the sensitive variables 0 and u, they have same 
distribution patterns. Meanwhile, with respect to the insensitive v, one is the others’ 
mirror image. 

Because G* ( ^ ( 0 , u, v) and G ^ , (0, u, v) are three-dimensional, it is impossible 
to plot their profile in one figure. We calculate several two-dimensional distributions of 
G a( Ct)) ( a > /*o. 7) over the Euler angles a and 7 at 296 K obtained with the fixed 0 O = 5, 22, 
38, and 50 degrees, respectively, and present their corresponding three-dimensional plots in 
Fig. 1. From the figure, one can easily see that the magnitudes of G * ( ( - T)) decreases very 
fast as 0 O increases. More specifically, for 0 O = 22, 38, and 50 degrees the magnitudes 
decrease by about one order each. In addition, these G* ( exhibit symmetry with respect 
to the axes (a + j)/2 and (a — y)/2. In order to show the profile of G* ( ^ , at 296 K over 
the two sensitive variables 0 and u, we calculate its average over v and present the 
resulting two-dimensional distributions in Fig. 2. One has to imagine that the profile 
shown in Fig. 2 extends along another dimension, i.e., the v axis which is perpendicular to 
the 0 — u plane and is missing in the figure. Furthermore, the range of extension along the v 
axis varies from the minimum 0 at u = 0 and 2ir to the maximum 27r at u = ir. As shown 
in Fig. 2, there are five sharp peaks located along the u axis at u = 0, tt/2, 7 r, 3 tt/ 2, and 27 t, 
respectively, and they are symmetric with respect to the plane u = n. The magnitudes of 
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these peaks decrease very fast as 0 increases. We note that in contrast with Fig. 1, a 
logarithmic coordinate is used to plot the magnitudes. From Figs. 1 and 2, one can 
conclude that the peak at u = 7r is dominant. We do not present the profiles of G * ( ^ } for 
other temperatures, but simply mention that they have similar patterns, but the peaks 
become lower and wider as the temperature decreases. On the other hand, it is 
unnecessary to present similar figures for G'^, because one can easily obtain them from 
Figs. 1 and 2. In fact, Figs. 1 and 2 are also applicable for G^ ( except that one has to 
switch the labels a and 7 in Fig. 2. Finally, we note that because G * ( ^ (0, u, v) and 
G * ( (0, u, v) are independent of the potential, it is wise to calculate them first and store 

them in files. Then, when one carries out repeated calculations for and X.(u) to 
optimize the potential models, one does not need to evaluate the values of G a(^)(^> u > v ) 
and G^ ( ^ (0, u, v) again. In addition, using the interpolation method one can easily 
obtain their values for a random selection of 0, u, and v from these input files. Otherwise, 
one has to independently evaluate them about 10 7 times in the Monte Carlo calculations. 

D. Contributions from Frequency Detuning 

So far, our discussion has been focused on how to apply the coordinate 
representation for calculating / + (w) and In order to calculate <T + (cj) and X.(u), one 

has to go further by adding contributions from the frequency detuning. We briefly explain 
a method used to obtain the second term of Eq. (22) associated with X + (u>). In comparison 
with evaluating the first term <T + (u>), the only additional obstacle is that the integrand 
contains a derivative x’ . (<*>). In practice, except for special cases in which analytic 

k+>n 

expressions for this derivative are available, to evaluate values of 7^ (w) is much more 

difficult than 7 . (w) because a numerical subroutine is required to obtain the former from 

k+>n 

the latter. Given the fact that there are about 10 7 random selections in the Monte Carlo 
calculations, this means that this subroutine must be called 10 7 times. Fortunately, except 
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for this part, the other parts of the second term do not depend on the frequency. 

Therefore, it is better to reverse the order of the derivative and the summation (i.e., 
integral) operations. In other words, instead of jjj- ^ ^ [ui — 5(0^* 4- o> n )] y’ . (u/), we can 

k+ n 

calculate a new term given by hll [w — 3(0^+ + w n )] x k . (w)- The results obtained are 

n 

a function of the frequency represented by a set of values of the integrations and a set of 
corresponding frequencies. Then, with the numerical subroutine, one is able to obtain the 
derivatives which are just the second term we want to calculate. As explained above, the 
subroutine is called only once. With this technique, the costs to calculate the second term 
of Eq. (22) is comparable to the first. 

The same method used for X 4 (u) explained in Sec. II. C is also applicable for 
evaluating this new term. We do not repeat a detailed discussion, but simply mention 

T ♦ 

things which are different. In this case, a new set of coefficients B , can be introduced 

A I X 

whose expression is almost the same as , given by Eq. (39) except a factor of {u; — 

An 

^[E(j^Tj) — EQjtJ + E(j^r£) — E(j 2 r 2 )]} is inserted into the inside of summation loops over 
the indices j t , t v j 2 , r 2 , j’ t , r\, j 2 , and r 2 . Similarly, one can introduce a set of coefficients 
, associated with the second term of Eq. (24) for X_(v) whose expression is the same as 

lx lx 

except the summation over jj and t \ is limited to a range with EQ’.r’.) < EQ.r,). We 

JEx lx 

L+ l “ L + r “ T + 

note that in comparison with A KR , and A KK „ B KK , and B KK , have similar symmetries B KR , 

T ♦ T- T “ T+ T " 

= B_ k _ k , and B kk , = B_ K _ K „ except B KK , = — B K , K< Then, one can introduce two 

T ♦ T - 

three-dimensional distributions associated with B KK , and B KK „ respectively, and store 
their values in two input files the same way as G* ( ^ (0, u, v) and G^ ( (0, u, v). 

Finally, by comparing these two distributions, their patterns are closely related each other 
as shown by Eq. (41) except one has to add a minus sign on the left side since their values 
become opposite. 

As expected, to evaluate B , and B , requires considerable CPU time, and the 

AA Aa 

T ♦ T - 

costs are approximately the same as for A , and A ,. In fact, for specified L , K and K’, 

KK KK 


20 


one does not need to calculate the four coefficients separately and one can evaluate all of 
them simultaneously. 

E. A Monte Carlo Calculation of 11— Dimensional Integrations 

As an example, we explain how to calculate l’+(w) from Eq. (23) in detail. In the 
coordinate representation, the summation of over ( and r/ becomes a 

11-dimensional integration of H^ T) (w)G ( + ( - T)) over the Euler angles /3 a <-, 7 a <-, a b(; , 0 bC> , 7 b( -, 
a aTi> 7an» a bTi> ^bn> and 7br| in which the first five (including a aC = 0) specify the initial 
orientations of the system and the last six specify the final ones. We note that due to the 
rotational symmetry of the whole system, one can always assume a al » = 0. For such high 
dimensionality, the Monte Carlo method is the only way to evaluate the integrals. 

It is well known that in the Monte Carlo computation, it is important to distinguish 
the sensitive and insensitive variables of the integrand, and to incorporate this into the 
integration variables since this enables one to tailor the important sampling and to reduce 
the variance dramatically. In the present case, the integrand is a product of and 

G(*;ti) • With respect to their variables, the former is a smooth function as shown by Eq. 
(9), but the latter’s values vary wildly and could differ from each other by many orders of 
magnitude. This means that the sensitivity of integrand is mainly determined by G^, , 
or more specifically, by G a( and G b( . We note that G a( and G b( ^ , depend on 
relative orientations between the initial and final positions of the absorber and disturber 
molecules, respectively. Therefore, it is proper to represent the final orientations of the 
system labeled by tj in terms of the body— fixed frames instead of the space— fixed frame. 

The body— fixed frames introduced here are those attached to the two molecules at their 
initial orientational positions. For the asymmetric— top molecule, one chooses a 0^, 7^, 
a (CT)) , , and 7^) instead of 0^, 0^, 7^, o^, 0^, and 7^ as variables. However, similar 

to the behavior of G b( ^ explained in our previous study, 3 the sensitivity of G a( ^ } on its 
two variables and 7^) is interwoven such that neither is a sensitive or insensitive 
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variable, but their combinations u ( ^ T)) , and v ( ^ are. Therefore, in order to well 
characterize the sensitivity, a replacement of a { ( - T)) , ^ , and 7 ( ^ by ( - T)) , u ( ^ T)) , and 

v ( as variables in G * ( ^ is appropriate. The same conclusion is also true for G b( ^ . 
This is a further step necessary to evaluate 11— dimensional integrations because, in 
comparison with 7 or 9— dimensional ones, not only the dimensionality becomes higher, but 
also the distributions of the integrand become more nonuniform. However, in order to 
incorporate these new choices for the integration variables, one has to pay extra attention 
to their ranges. More specifically, since u ( ^ varies from 0 to 27 t and v ( ^ varies from — 
u ( to u { ^ T)) when 0 < u < ir and from u ( ^ — 2ir to 2x - u ( ^ T)) when 7 r < u < 27 r, the 
integration volume becomes a lozenge-shaped area. Because the algorithm VEGAS 10 is 
designed for carrying out integrations over rectangular volumes, one cannot incorporate the 
integration variables directly. Therefore, we have modified VEGAS such that the new 
version enables one to evaluate integrations over a volume containing a lozenge— shape 
area. Then, with respect to the sensitivity, the dependence of integrand on all integration 
variables is well characterized. This enables one to fully exploit the power of the Monte 
Carlo method. As a result, numerical tests show that with a few more random selections 
than before, one is able to evaluate the 11-dimensional integrations successfully. 

The above discussion is also applicable for evaluating other 11— dimensional 
integrations, i.e., those for and for the frequency detuning correction terms of X+(u) 
and because their integrands have similar features. Numerical tests show that with 
the modified version of VEGAS, we are also able to evaluate these as well. 

F. Potential Models 

Based on the progress mentioned so far, we are able to calculate converged line 
shapes for H 2 0 — H 2 0 without or with the frequency detuning correction, i.e., <T + (u/) and 

or /? + (w) and <T.(w), from potential models provided unless they are too complicated 
(e.g. those consisting of several decades terms and parameters). There are several potential 
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models available in literatures, such as HF 11 , Cl 12 , Watts 13 , RWKl 14 , and RWK2 model 14 . 
We have tested some of them to calculate the line shapes and the corresponding absorption 
coefficients. Unfortunately, it turns out that the results obtained from these models 
predict too much absorptions in comparison with experimental data. It has been shown in 
our previous study 2 ’ 3 that the far— wing shape is very sensitive to the angular gradients of 
the potential. The reason is that contributions to the line shape come mainly from energy 
contour areas at which large angular gradients of the potential are exhibited while the 
potential values themselves are relatively small or even negative. We expect that such 
features would not fully manifest their effects on other physical measurements on which 
these models were developed. Therefore, the failure of these models is not surprising. 

As an alternative, we assume that the potential for H 2 0 - H 2 0 consists of a 
Coulomb interaction represented by a site-site model, a short range repulsive interaction 
represented by another site— site model, and an isotropic attractive dispersion interaction 
proportional to — 1/r 6 . For each H 2 0, there are two positive point charges + q located at 
the H atoms and one negative point charge - 2q at a position along its symmetry axis a 
distance d from the 0 atom. (A positive value of d means the charge — 2q is located on the 
same side as the center of mass.) In addition, we assume that there are three repulsive 
force centers: two located at the two H atoms and one at the 0 atom, and the repulsive 
interactions have a form Ay e r ij/^ij, where ry are distances between force centers, and 
Ay and py are adjustable parameters. Accordingly, the potential V(r, fl a , fl b ) considered 
here is given by 

V(r, 0., 0„) = l l 3191 + l l Ay r - f, , (42) 

ieajeb 1J i€aj€b 

and it contains 7 adjustable parameters: A Q0 , p QQ , A 0}{ , p 0JJ , A HH , p HH , and B. The values 
of q = 0.60 (a.u.) and d = 0.4991 (a.u.) are determined such that they match the well 
known dipole moment value 1.8546 D and yield quadrupole moments $ cc = — 2.319 D A, 

0aa = 2.635 D A, and 0 bb = — 0.316 D A, which are reasonable in comparison with 
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experimental values 0 CC = - 2.50 D A, 0 aa = 2.63 D A, and 0 bb = - 0.13 D A . 15 We adopt 
the potential form given by (42) and search for a set of potential parameters that enables 
us to obtain satisfactory results for several properties of a dilute H 2 0 gas; e g. the 
absorption coefficient, the second virial coefficient, and the differential cross section. 

G. General Features and Numerical Results for the Line Shapes 

We have presented the line shape formulas applicable for both the pure rotational 
band and vibrational bands. But, in the present study, the numerical calculations are 
carried out for the former because not only is this the strongest band of H 2 0, but also most 
of the continuum absorption measurements are performed in its high-frequency wing. 

It is worthwhile to report general features of line shapes found from numerous test 
calculations before we present more specified results. First of all, we find that differences 
between X^ui) and X_(u) calculated from the same potential are always less than n um erical 
errors. The formulas used to get X t (u>) and X.{u) are the same, but the input files 
representing the two distributions G a( ^ (0, u, v) and G a( ^ (0, u, v), respectively, are 
different. However, these two distributions differ from each other only slightly. In fact, as 
explained above, they exhibit the same profiles over two sensitive variables 0 and u, and 
are mirror images over the insensitive v. In cases that the potential contains cyclic 
coordinates, it is easy to show analytically that these two distributions must yield the same 
line shapes. For more general cases, given the fact that G a( (0, u, v) and G a( ^ , (0, u, 
v) differ from each other as discussed above, we suspect that effects resulting from these 
differences could cancel out in the averaging processes when the integrations are evaluated 
by about 10 7 random selections. In the present study, we do not seek a general and 
rigorous proof of this finding, rather we assume it. Then, we can draw the important 
conclusion that by not considering the frequency detuning correction, the line shape 
obtained from the formalism satisfying the detailed balance is symmetric. In other words, 
if one does not consider the frequency detuning correction, one only needs to introduce one 
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line shape function. 

Second, we find that the calculated contributions from the frequency detuning 
associated with Xj^d) are opposite to those associated with X.(d), and the differences 
between their magnitudes are always less than numerical errors. This indicates that these 
correction contributions have the same magnitudes, but with different signs; this is not 
surprising since we have already noted our finding for X+(d) and X.(d). In general, the 
values of the correction contribution for XJ^d) are positive and those for X.(oj) are negative 
for u > 0 and vice versa for u> < 0. In addition, we find that the correction contributions 
are significant, especially at high frequencies, but they tend to become negligible near the 
line center. 

Now, we are ready to present some results. We find that by adopting A Q0 /k = 1.05 
» 10’ K, p 00 = 0.245 A, A 0K /k = 2.0 < 10> K, p m = 0.36 A, A >|( = 4.0 . 10 2 K, p m = 

0.46 A, and B/k = 9.0 * 10 3 K, one is able to obtain X,(u) and X.(u>) such that the 
calculated absorption at 296 K can fit the experimental results in the spectral region 300 - 
1100 cm" 1 very well. We note that because this window region is located at the 
high-frequency side of the pure rotational band, the calculated absorption are mainly from 
the contributions of X+(u). The two line shapes X+(d) and X.(u ) at 296 K, together with 
their mean X(u), are plotted in Fig. 3. As shown in the figure, the magnitudes of X + (d) are 
larger than X.(d) and the relative gaps between them increase significantly as u increases. 
The corresponding self-broadened absorption coefficients in the spectral region 300 - 1100 
cm" 1 based on HITRAN 92 data are plotted in Fig. 4, together with the experimental 
results of Burch et al. 16 and some recent measurements of Cormier et al. 17 Using this 
potential, we calculate the second virial coefficients at several temperatures 18 and compare 
them with experimental data 19 in Fig. 5. In addition, the calculated differential cross 
sections together with experimental data 20 are plotted in Fig. 6. We note that all these 
physical quantities are associated with dilute water vapor, and as shown by these figures, 
the agreements between the theoretical predictions and the laboratory measurements are 
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good. 

Based on the same potential model, we calculated the two line shapes X.(u) and 
X{u) for several temperatures ranging from 220 K to 330 K which are of interest in the 
atmospheric applications and the corresponding absorption coefficients. Some of X+(u) 
obtained for 220, 240, 260, 280, 300, and 320 K with frequencies raging from — 1600 cm* 1 to 
1600 cm* 1 are presented in Fig. 7 and the corresponding X,(u) are their reflections about 
the u = 0 axis. As shown in the figure, these X.(u) and X.(u) are asymmetric and their 
magnitudes increase as the temperature decreases. We note that these line shapes X+(u) 
and X.(u) presented here do not include a factor 1/u 2 as shown in Eq. (13). If one wants to 
compare them to other line shapes (e.g. a Lorentzian), one has to multiply them by the 
factor 1/u 2 . Also, the magnitude of X^v) or X.(u) as w-» 0 should approach the Lorentzian 
halfwidth, although the present theory is not valid in this limit. It is clear, however, that 
the theoretical shapes are super— Lorentzian for the displacements up to around 400 cm' 1 , 
and then become sub— Lorentzian for larger displacements. In fact, the line shape must 
approach zero at least as fast as an exponential. This can be shown from the analyticity of 
the correlation function; the successive derivatives of C(t) in the t = 0 limit correspond to 
moments of the line shape in frequency space. Because the derivative are all finite, this 
implies that all the moments of the line shape are also finite and, therefore the line shape 
must approach zero faster than any inverse power of u. In Fig. 8, we present all the 
calculated absorption coefficients in the window region 600 — 1250 cm' 1 . We note that for 
a specified frequency u, we exclude completely any contribution from lines that are within 
[u — 25 cm* 1 , u + 25 cm* 1 ] in the calculations. As shown by Fig. 8, the negative 
temperature dependence of self-continuum is clearly demonstrated. 

III. DISCUSSIONS AND CONCLUSIONS 

In comparison with our previous studies on the far— wing line shape, there are 
several important advances which have been made in the present study. First of all, by 


/ 


26 


dearly distinguishing the sensitive and insensitive variables and by modifying the Monte 
Carlo subroutine used previously to handle integrations whose volume is not rectangular, 
the effectiveness of the important sampling is enhanced significantly. As a result, by 
accounting for random selections the order of 10 7 (which is comparable to or slightly more 
than that required for lower dimensionality cases), one is able to evaluate 11— dimensional 
integrations. 

Second, we have carried out numerical calculations based on the formalism which 
satisfies the detailed balance principle exactly and has adiigher accuracy in the short-time 
limit. As expected, this increases the difficulty because the evaluation of the integrands 
requires more calculations. As shown by Eqs. (3) and (4), within this formalism the dipole 
moment operator appears in formulas in such a way that it is always sandwiched by the 
density matrix. However, these two operators have different characters: one depends on 
the coordinates only while the other contains the differential operators. No matter what 
kind of representation is chosen, the sandwiched operators require more loops to evaluate 
their values. As in our previous studies, we use the coordinate representation because it 
enable us to include more populated states; we also introduce the distribution functions and 
store them in files because it enable us to obtain values of the integrand with less CPU 
time. But, to calculate these distributions with a high cut-off j max requires long CUP 
times. We note that for temperatures of interest in atmospheric applications, j,^ = 23 is 
enough, but for higher temperatures, a higher j max would be necessary. With these input 
files, to obtain values of the integrand becomes relatively easy. Combing the technique 
mentioned above to reduce the number of random selection to the order of 10 7 , one is able 
to complete the evaluation of 11— dimensional integrations and to obtain a line shape in one 
day with one workstation. 

Finally, by carrying out band averages in a more sophisticated way, we can 
calculate the two line shape functions. In addition, the effects resulting from the frequency 
detuning have been taken into account in the averaging processes. We note that except for 
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the simplest system, e.g. C0 2 — Ar, 8 the latter has not been done previously. In fact, most 
of calculations are similar to those without including the frequency detuning correction, 
except the integrand contains a derivative. As explained above, we have developed a 
technique which enables one to overcome this obstacle and calculate the frequency detuning 
corrections. 

Based on numerous test calculations of line shapes during the course of study, we 
can draw several important conclusions. It is necessary to consider the frequency detuning 
correction because the effect on the line shape is significant, especially in the 
high-frequency region. However, this effect shows up only when one distinguishes the 
positive and the negative resonance lines and carries out the band average over them 
separately. In other words, we have demonstrated that the line shapes obtained from the 
two band averages is asymmetric and we have found that this asymmetry results from the 
frequency detuning effect, or more specifically, from the distribution of lines within the 
band of interest. In addition, we find that for the pure rotational band the magnitudes of 
<T+(w) are significantly larger than X.(u>) for u > 0 and vice versa for u < 0, and these gaps 
become larger as w increases. We note that the conclusion concerning the asymmetry 
claimed here is applicable for the band average line shapes and has nothing to do with 
individual lines. 

Finally, we would like to make a few of comments on the vibrational bands. We 
expect the main conclusions about the line shapes for the pure rotational band would 
remain true, but some different features could show up because the contributions to X+(u) 
and X.(u) from the frequency detuning terms depend strongly on the band structure. We 
expect that the more unevenly and the more widely the lines are distributed within the 
band of interest, the contributions come from the frequency detuning corrections will 
increase and the more / + (w) differs from X.(u). It is well known that the line distributions 
of the vibrational bands are quite different from that for the pure rotational one. The 
former’s lines are, more or less evenly, located on both sides of the band centers, but the 
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latter’s are always on the high-frequency side because the band center is zero. We note 
that according to the definition, the average positive resonance frequency u> introduced here 
is not the band center. But these two are very close for the vibrational bands and are quite 
different for the pure rotational band. We expect that for the vibrational band the 
frequency detuning corrections could become smaller. This means that XJ^u) and X_(u) 
differ from each other by smaller amounts. Meanwhile, we could not draw any conclusion 
about which magnitude is larger because this is related to the special structure of the pure 
rotational band. It has been known for years that the line shapes for C0 2 derived from 
experimental data are asymmetric and they are not the same for different bands. 21 ’ 22 So 
far, there has been no theoretical explanation why the different bands have different shapes 
without assuming that the interaction depends sensitively on the vibrational quantum 
numbers which seems unlikely to be true. We think that both from the theoretical and 
practical points of view, to investigate the lines shapes for different bands is an interesting 
subject to pursue. 
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LIST OF FIGURES 

Fig. 1. The two-dimensional distribution of G ^ (a, 0, 7) of H 2 0 over the Euler angles 
a and 7 obtained at T = 296 K for j,^ = 26. The values of the Euler angle 0 is 
fixed and the four plots presented here correspond to 0 = 5, 22, 38, and 50 degrees, 
respectively. 

Fig. 2. The two-dimensional distribution of G * ( ^ Tj) (0, u, v) of H 2 0 over the two sensitive 
variables 0 and u obtained at T = 296 K for = 26. This distribution results 
from the averaging G * ( ^ T)) (0, u, v) over the one insensitive variable v. In contrast 
with Fig. 1, a logarithmic coordinate is used for the G* ( ^ , (0, u, v) axis. 

Fig. 3. The self— broadened far— wing line shape of H 2 0 (in units of cm' 1 atm' 1 ) as a 
function of frequency ui (in units of cm' 1 ) for T = 296 K. The dashed curve 
represents <T + (o;) calculated from the positive resonance line average and the dotted 
curve represents X.( u ) calculated from the negative resonance line average. The 
solid curve is X(u) which is the mean of <T + (u>) and 

Fig. 4. The calculated self— broadened absorption coefficient (in units of cm 2 molecule' 1 

atm' 1 ) at T = 296 K in the 300 — 1100 cm' 1 spectral region is represented by a. For 
comparison, the experimental values of Burch et. al are denoted by + and those 
from Cormier et al. are denoted by □. 

Fig. 5. The calculated second virial coefficients (in units of cm 3 mol' 1 ) as a function of 

temperature is represented by a solid line. The experimental data are denoted by □. 

Fig. 6. The calculated differential cross section (in arbitrary units) as a function of the 

laboratory scattering angle 0 is represented by a solid line. The experimental values 
of Duquette are denoted by □. 

Fig. 7, The self— broadened far— wing line shapes <T + (w) of H 2 0 (in units of cm' 1 atm' 1 ) as a 
function of frequency u (in units of cm' 1 ) obtained for T = 220, 240, 260, 280, 300, 
and 320 K; these are represented by 6 curves in order from top to bottom. The 
frequency u varies from — 1600 cm' 1 to 1600 cm' 1 . The corresponding X.(u) are 
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reflections of X+(u>) about the u = 0 axis. 

Fig. 8. The self— broadened absorption coefficient (in units of cm 2 molecule* 1 atm* 1 ) in the 
window region 600 - 1250 cm* 1 calculated for T = 220, 230, 240, 250, 260, 270, 280, 
290, 300, 310, 320, and 330 K in order from top to bottom. A cut-off 25cm* 1 is used 
to exclude completely any contribution of lines that are closer than this limit. 
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